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We study phonons in a dynamical chain of ions confined by a trap with a time-dependent (axial) 
potential strength and demonstrate that they behave in the same way as quantum fields in an 
expanding/contracting universe. Based on this analogy, we present a scheme for the detection of 
the analogue of cosmological particle creation which should be feasible with present-day technology. 
In order to test the quantum nature of the particle creation mechanism and to distinguish it from 
classical effects such as heating, we propose to measure the two-phonon amplitude via the 2"^* red 
side-band and to compare it with the one-phonon amplitude (1^' red side-band). 
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Introduction The theory of quantum fields in curved 
space-times (see, e.g., [Ij) comprises many fascinating 
and striking phenomena - one of them being the creation 
of real particles out of the (virtual) quantum vacuum fluc- 
tuations by a gravitational field. These effects include 
Hawking radiation given off by black holes as well as cos- 
mological particle creation. A very similar mechanism - 
the amplification of quantum vacuum fiuctuations due to 
the rapid expansion of the very early universe - is (ac- 
cording to our standard model of cosmology) responsible 
for the generation of the seeds for cosmic structure for- 
mation. Hence, even though these effects are far removed 
from every-day experience, they are very important for 
the past and the future fate of our universe. 

Therefore, it would be desirable to render these phe- 
nomena accessible to an experimental verification. Prob- 
ably the most promising way for achieving this goal is to 
construct a suitable analogue which reproduces the rele- 
vant features (such as the Hamiltonian) of quantum fields 
in curved space-times. Along this line of reasoning, pro- 
posals based on the analogy between phonons in dynam- 
ical Bose-Einstein condensates and quantized fields in an 
expanding/contracting universe have been suggested 
Unfortunately, the detection of the created phonons in 
these systems is rather difficult (see, however, [1]). 

On the other hand, the detection of single phonons in 
ion traps via optical techniques is already state of the art 
in current technology - which suggests the study of this 
set-up instead. In this Letter, we shall derive the analogy 
between phonons in an axially time-dependent ion trap 
and quantum fields in an expanding/contracting universe 
and propose a corresponding detection scheme for the 
analogue of cosmological particle creation. A similar idea 
has already been pursued in but the proposal pre- 
sented therein goes along with several problems, which 
will be discussed below l5l| . 

Cosmological particle creation Let us start by briefiy 
reviewing the basic mechanism of particle creation in an 
expanding/contracting universe. For simplicity, we con- 



sider a massless scalar field (j) described by the action, 
see, e.g., [ij {h — c=\ throughout) 

A^\j d^x^\ [(9;.<^)fl^'^ (9.0) - C^<i?] , (1) 

where q^^" denotes the metric and g = detjg^j,} its de- 
terminant. Furthermore, a scalar field can be coupled 
to the Ricci (curvature) scalar 9^ via a dimensionless pa- 
rameter C (e.g., conformal coupling C = 1/6, cf. [H). A 
spatially flat universe can be described in terms of the 
Friedman-Robertson- Walker metric 

ds^ = a'^{t)dt^ ~a'{t)dr^ , (2) 

with the time-dependent scale parameter a{t) corre- 
sponding to the cosmic expansion/contraction. Here we 
have chosen a slightly unusual time-coordinate t related 
to the proper time r via dr — a^{t)dt in order to simplify 
the subsequent formulae. After a normal-mode expan- 
sion, the wave equation reads 

4>k + [a\t)k^ + Ca^(i)«(0] (bk = 0, (3) 

i.e., each mode k just represents a harmonic oscillator 
with a time-dependent potential a^{t)k'^ + (^a^{t)yi{t). As 
long as this external time-dependence of the potential is 
much slower than the internal frequency of the oscilla- 
tor, the quantum state will stay near the ground state 
due to the adiabatic theorem. However, if the external 
time-dependence is fast enough (i.e., non-adiabatic) , the 
evolution will transform the initial ground state into an 
excited state in general - which is the basic mechanism 
of cosmological particle creation. In this case, the initial 
vacuum state |0) = \ip(t | — oo)) containing no particles 
Vfe flfc |0) =0 evolves into a squeezed state 

IV(noo)) = exp|x^6 4atfe-h.c.| |0) 

= \o)+j2^k\ik,i-k)+oiek), (4) 
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which does contain pairs of particles |lfc,l-fe)- The 
squeezing parameter for each mode is determined by 
the solution of Eq. ^ and thus by the time-dependence 
of a'*(i)fc^ as well as (a^(t)^R{t) and governs the number 
of created particles per mode 

{hk) = {^p{t T alak mt T (»)) = sinh^defel) . (5) 

Ion-trap analogue Assuming a strong radial confine- 
ment of the ions, we consider their axial motion only. In 
a time-dependent harmonic axial potential described by 
the oscillator frequency tOaxit), the position qi of the i-th 
ion obeys the equation of motion 
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sign(z - j) 



(6) 



where the factor 7 encodes the strength of the Coulomb 
repulsion between the ions. Assuming a static situation 
initially, the classical solution to the above equation can 
be obtained via the scaling ansatz qi{t) = b{t)q^, where 
q^ are the initial static equilibrium positions, leading to 
the evolution equation for the scale parameter b{t) 
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(7) 



In order to treat the quantum fluctuations of the ions 
(leading to the quantized phonon modes), let us split the 
full position operator qi{t) for each ion into its classical 
trajectory b(t)qf and quantum fluctuations 5qi{t) 



(8) 



Since these fluctuations dqi{t) are very small for heavy 
ions, we may linearize the full equation of motion ^ 



+ ^Lit) 5q^ = 



1 



bHt) 



^M^jSqj 



(9) 



with a time-independent matrix My arising from the 
Coulomb term in fH]). Diagonalization of this matrix 
(normal-mode expansion) yields the phonon modes 



b^t) 



Sqi^ = 0, 



(10) 



labeled by k. The time-independent eigenvalues ui'^ > 
of the matrix Mij determine the phonon frequencies. The 
lowest mode is the center-of-mass mode corresponding 
to a simultaneous (rigid) motion of the ions. Since the 
ion distances are fixed, the Coulomb term does not con- 
tribute in this situation = 0. The next mode is the 
breathing mode with ujf. — 2^2^(0). Comparing Eqs. ([7]) 
and (fTO]) . we see that this mode exactly corresponds to 
the scaling ansatz, i.e., the ion cloud expands/contracts 
linearly. Hence this is the only mode which can be ex- 
cited classically (for a purely harmonic potential). I.e., 



without imperfections such as heating, phonons in the 
other modes can only be created by quantum effects. 

Comparing Eqs. ([3]) and (jlOp and identifying 0^ with 
Sq^, we observe a strong similarity: The wavenumber 
in ^ directly corresponds to cj2 in ([TU|) and the scale fac- 
tors a(t) and b{t) enter in a similar way. However, an ex- 
panding universe is analogous to a contracting ion cloud 
and vice versa. In the mode-independent terms, the axial 
trap frequency t^ax acts like the Ricci scalar fH. Interest- 
ingly, both are related to the second time-derivatives of 
the corresponding scale factors. 

In view of the formal equivalence of Eqs. ([7]) and (|10p . 
we obtain the same effects as in cosmology - in particular, 
the mixing of creation and annihilation operators 



(11) 



which can be expressed in terms of the Bogolubov coef- 
ficients satisfying \a'^\ — |/32| = 1. Note that the above 
relation is just the operator representation of the squeez- 
ing transformation in Eq. (jH) with \(3^\ sinh(|^/j|). 

Detection scheme In the following we describe how 
to realize the proposed experiment by applying oper- 
ations closely related to those implemented on qubit- 
ions in quantum information processing We focus 
on initializing the system, simulating the non-adiabatic 
expansion of space, performing the read-out of the final 
state (particle- or phonon-number distribution) and dis- 
tinguishing it from a classically describable outcome, for 
example caused by thermal heating. To perform a first re- 
alization, we will confine one single earth alkaline atomic 
ion to the axis of a linear radio- frequency trap [l3 | similar 
to that described in 0, Q ■ The required simulation basis 
can be composed by a ^Si/2 electronic ground state level 
of 25Mg+, here the state |F = 3;m/ = 3) = and the 
associated harmonic oscillator levels \n) related to the 
axial harmonic confinement, as depicted in Fig. [TJ At 
the start of each experiment, the ion will be laser cooled 
close to the ground state of the axial (external) motion 
and optically pumped into the electronic (internal) state 
li) @]- Then we will decrease adiabatically the axial con- 
finement and subsequently reset it non-adiabatically to 
its initial value (as already proposed in [l3| in another 
context). Since the ground state wave function of the 
ion cannot adapt to the restored stiff confinement (non- 
adiabatic case) , it will oscillate symmetrically around the 
minimum of the final trapping potential [ll[, i.e., without 
populating odd motional states. As shown above, this 
non-classical oscillation is to be described via a squeezed 
state (see also 12]) depicted in Fig. [2] 

We propose, in addition to established schemes de- 
scribed in [l2| or [l3|, for example, an alternative method 
to distinguish classical noise (such as the initial ther- 
mal distribution or heating during the process) from a 
squeezed state generated by quantum effects considered 
here. In order to read out the final motional state, we 
will first couple it (via suitable lasers) to two internal 
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FIG. 1: Schematic of the relevant energy levels (not to scale) 
of one ^^Mg^ ion. Shown are the ground-state hyperfine levels 
supplying the two internal states (|J,) and ||)) and the first 
three equidistant harmonic oscillator levels |n), related to the 
harmonic axial confinement in a linear ion trap. Typically, the 
energy splitting of the motional levels and the Zeeman shift 
induced by an external magnetic field are of the same order 
of magnitude within 1-10 MHz, therefore much smaller than 
the Hyperfine splitting of 1.8 GHz, the fine structure splitting 
of 2750 GHz and the optical transition frequency of the order 
of lO^^Hz. We depict the resonant transition state sensitive 
detection named (d) and the relevant types of ofi^- resonant 
(~ 100 GHz) two-photon stimulated Raman transitions (a,b 
and c) described in the text. 

states of the ion. Besides the electronic ground state 
\F = 3; TO/ = 3) = II), the second internal state to 
span a two-level system (analogous to a qubit) is im- 
plemented via a second hyperfine state of ^^Mg+, \F ~ 
2, TO/ = 2) = It), separated from the state ||) by the 
hyperfine splitting Wo — 27r x 1.8 GHz. We will accom- 
plish the coupling of the two internal states ||), |t) and 
the motional states \n) via two-photon stimulated Ra- 
man transitions Q requiring two laser beams (A « 280 
nm), with wave vector difference Afe = k-^ — aligned 
along the trap axis z (|AA:| = \/2 x 27r/A = 27r/Aefi). Via 
detuning the frequency difference uj2 — loi from the hyper- 
fine splitting ojo ± 2'Kmvz by integer multiples to of the 
axial trapping frequency Vz, we may drive the carrier- 
transition (to = 0) or the first- (to = 1) and second- 
(to — 2) sideband transitions respectively. Note that 
the spectral resolution of the two-photon stimulated Ra- 
man transition is independent on the natural line width 
r = 27rx43 MHz of the resonant transition - but propor- 
tional to the inverse of the Rabi-frequency, adjusted via 
the intensities or detuning of the laser beams allowing 
for the resolution of the individual motional states sep- 
arated hy Vz -^T. In order to measure the population 
of the motional state \n = 2), we will drive a sequence of 
transitions (cf. Fig. [1]) , synthesized by a second-sideband 
(a) transition (||,n = 2)^ |'f,?T. = 0)) followed by a car- 
rier (c) transition (||,n = 0) |J,,n = 0)). The final 
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FIG. 2: Population of the motional states: On the left hand 
side, we show the experimentally realized thermal spectrum 
after cooling the system close to the ground state ((n)=0.05). 
On the right hand side, we show the simulated spectrum of 
the ((n)=l) squeezed state, assuming the starting distribution 
from the left. The almost complete absence of population in 
odd motional states emphasizes the non-classical character of 
this state and guides towards a simple possibility to distin- 
guish it from thermal and disturbed ones via measuring a 
much higher population in state \n = 2) than \n — 1). 

read-out (d) described below is internal-state dependant 
and provides us with the population in state ||). After 
the sequence (a,c,d) of transitions this is almost exclu- 
sively equivalent to the population of the motional state 
|n = 2) (because the probability of even higher excita- 
tions |n > 3) is expected to be much smaller and their 
Rabi frequencies are also different). This result can then 
be compared with the outcome after a first red-sideband 
(b) transition (|J,,n=l)^|t,n = 0)) followed by a car- 
rier (c) transition (||, n = 0) ||, n = 0)), providing the 
probability of motional excitation \n = 1). As soon as 
we deduce a higher probability for motional state |n = 2) 
than for motional state \n = 1), we show the incompat- 
ibility with classical effects such as a thermal distribu- 
tion and get strong evidence for the non-classical effect 
of squeezing [Tl ]. 

Finally, we have to read out the internal state effi- 
ciently. To this end, we apply an additional resonant 
laser beam (d) , tuned to a cycling transition , coupling 
only state ||) resonantly to the P3/2 level and providing 
spontaneous emission at rates of > 10 MHz. This allows 
to distinguish the "bright" ||) from the "dark" |t) state 
with high accuracy, even at a low detection efficiency 
(due to the restricted solid angle etc.). 

Envisioned results The above mentioned sate of the 
art techniques allow to cool the axial motion close to 
the ground state (n) « 0.05 [3] (see also d 0]) and 
to optically pump into the down state ||) with 99% Q 
or even higher fidelity [^, 0]. First experiments show 
a possible non-adiabatic variation of the axial motional 
frequency i^z between 200 kHz and > 2 MHz with a re- 
lated rise time of the order of one micro-second, which 
is sufficiently fast compared to the oscillation period of 
the lower frequency. Numerical simulations (based on 
measured temporal variation curves) indicate that we 
should be able to transfer approximately 20% of the mo- 
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tional state population from the ground state |n = 0) 
into state |n = 2), which corresponds to a squeezed state 
with (h) « 1. Starting with a thermal distribution with 
(n) « 0.05 instead of the exact ground state \n — 0), 
there will also be a small final population (a few percent) 
of the state \n = 1), see Fig. [21 However, this resid- 
ual effect will be significantly smaller than the \n — 2) 
population such that the signatures of squeezing can be 
measured as described above. 

Conclusions Since the state of the art fidelities for the 
carrier and sideband transitions as well as the state sen- 
sitive detection exceed 99% [1, Q , the initialization and 
measurement of the system can be provided with high 
accuracy. In order to benefit from these operational fi- 
delities, the main task will be to minimize classical distur- 
bances. For example, we have to carefully balance the ap- 
plied voltages for confinement during their non-adiabatic 
changes to prevent classical excitation of the axial mo- 
tional mode. In comparison to some other experiments 
with ion traps, the requirements for the present proposal 
may be a bit easier to achieve because the duration of the 
experiment will be short (around 3 ms) compared to the 
inverse of the thermal heating rate for motional quanta 
inside the trap (<0.005 quanta/ms [3]) and because the 
thermal and the squeezed motional spectra show funda- 
mentally different characteristics, see Fig. [2j It should 
also be emphasized that it is impossible to resolve in- 
dividual motional states with pulse durations short com- 

gared to the inverse of their frequency difference (see also 
I and [Hi). This impossibility in resolution is related to 
the Heisenberg uncertainty principle AEAt > h/2 that 
allows to create pairs of particles (phonons) out of the 
vacuum (ground) state in first place. Increasing the sys- 
tem towards 8 modes (ions) might be feasible by this 
proposal with state of the art techniques [l6l |. further 
scaling might benefit from the technical progress driven 
by the attempts of the quantum information community. 

Apart from experimentally testing the analogue of cos- 
mic particle creation - which might ultimately allow the 
study of the impact of decoherence and interactions etc. 
- the investigation of non-adiabatic switching of trap- 
ping potentials and its influence on the quantum state 
on motion might also shed light on possible problems in 
schemes where the fast shuttling of ions in a multiplex 
trap architecture is required for scaling towards a univer- 
sal quantum computer. 
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